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Abstract 

It has been realized in recent years that coupled focusing lattices in accelerators and storage 
rings have significant advantages over conventional uncoupled focusing lattices, especially for high- 
intensity charged particle beams. A theoretical framework and associated tools for analyzing the 
spectral and structural stability properties of coupled lattices are formulated in this paper, based 
on the recently developed generalized Courant-Snyder theory for coupled lattices. It is shown 
that for periodic coupled lattices that are spectrally and structurally stable, the matrix envelope 
equation must admit matched solutions. Using the technique of normal form and pre-Iwasawa 
decomposition, a new method is developed to replace the (inefficient) shooting method for finding 
matched solutions for the matrix envelope equation. Stability properties of a continuously rotating 
quadrupole lattice are investigated. The Krein collision process for destabilization of the lattice is 
demonstrated. 

PACS numbers: 29.27.-a,52.20.Dq 
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I. INTRODUCTION 


The transverse focusing lattice is one of the few crucial subsystems in modern accelerators 
and storage rings. Most contemporary accelerators and storage rings are designed based on 


an uncoupled linear transverse lattice [l|, where the two degrees of freedom in the transverse 
directions are decoupled. Well-known analyses of the effects of weak coupling on stability 
properties have left the incorrect impression that the coupling between the .x-dynamics and 
^/-dynamics always results in instabilities or other deleterious effects. It has been realized 
recently that coupled lattices are not necessarily more unstable than uncoupled lattices. On 
the contrary, it is believed that coupled lattice can be more advantageous in comparison 
with conventional uncoupled lattices, especially for high-intensity charged particle beams 
(2 11]. This is because the parameter space for coupled lattices is much larger than that for 


uncoupled lattices, and one can explore the larger parameter space for a coupled lattice to 
optimize the lattice design. 

Of course, the most important consideration in lattice design is its stability properties. A 
thorough study of lattice stability requires one to distinguish two types of stability properties, 
spectral stability and structural stability 


12 


17]. The spectral stability of a linear periodic 


lattice is determined by the eigenvalues of the one-period map M of the lattice. If there exists 
a vector v such that M l v —> 00 as l —> 00 , then the map is spectrally unstable. Otherwise, it 
is spectrally stable. This is the most familiar stability property that is often analyzed. The 
structural stability of the lattice refers to the robustness of the spectral stability property of 
the lattice with respect to a structural perturbation, such as imperfections in the magnets, or 
misalignment of the beam-line. A lattice is structurally unstable if there exists a spectrally 
unstable lattice infinitesimally close-by. Otherwise, it is structurally stable. 

Unfortunately, our understanding of the stability properties of coupled lattices is far 
from comprehensive due to the lack of an effective mathematical tool to describe the coupled 
dynamics. For ID uncoupled dynamics, the de facto standard for parameterizing the focusing 


lattice is the Courant-Snyder (CS) theory |l|, which is mathematically elegant and directly 
linked to the physics of the beam dynamics. For coupled lattices, several parameterization 
schemes have been developed 181-426]. But none of these schemes is as effective for coupled 
lattices as the CS theory is for uncoupled lattices. Recently, we have developed a generalized 


Cour ant-Snyder theory for coupled lattices 


171. 127-34], which generalizes every important 
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aspect of the original CS theory to higher dimensions. Especially, the key components of 
the original CS theory, i.e., the envelope f un ction (or the j3 function) and the associated 
envelope equation are generalized into a matrix envelope function and the associated matrix 
envelope equation. 

In the present study, we apply the generalized CS theory to investigate the stability 
properties of coupled lattices. We prove an important proposition that a necessary condition 
for a periodic coupled lattice to be spectrally and structurally stable is that the generalized 
matrix envelope equation admits a matched solution. We also show how to apply the 


techniques of pre-Iwasawa decomposition 


35 


361] and normal form to construct a matched 


solution of the matrix envelope equation by simply solving the initial value problem once. 
This new method is of great value even for uncoupled lattices. Previously, one used the 
conventional shooting method to solve the initial value problem many times to search for 
a matched solution. Using the example of a continuously rotating quadrupole lattice, we 
illustrate in this paper how the lattice becomes spectrally unstable through an interesting 
process called the Krein collision. 

The paper is organized as follows. In Sec.QTl we describe the spectral and structural sta¬ 
bility properties of a generic linear periodic Hamiltonian system and the associated Krein 
collision. The generalized Courant-Snyder theory for coupled lattices is introduced in Sec. lIIIl 
and the connection between stability properties and matched lattice functions are discussed 
in SeclIVl The formalism developed is applied to study the stability properties of a contin¬ 
uously rotating quadrupole lattice in Sec.lVl 


II. SPECTRAL AND STRUCTURAL STABILITY PROPERTIES OF LINEAR 
HAMILTONIAN SYSTEMS 

The dynamics of a charged particle in a coupled or uncoupled periodic focusing system 
is completely specified by the one-period map. Because of the Hamiltonian nature of the 
dynamics, the one-period map is symplectic. For the linear focusing lattices considered in 
the present study, the one-period map is specified by a symplectic matrix M. In this section, 
we discuss the stability properties of M as a general symplectic matrix. The calculation and 
stability analysis of M for a specific choice of focusing lattice will be discussed in Secs. lIIIl 
and m using the generalized Courant-Snyder theory. 


3 










Let the dimension of M be 2n x 2 n, and let A/ (l = 1,2 n) be the eigenvalues of M. It 
is straightforward to prove that if A is an eigenvalue of a symplectic matrix, then its inverse 
1/A and its complex conjugate A are also eigenvalues. Then, the eigenvalue distribution can 
be divided into four categories: 


( 1 ) : All eigenvalues are distinct and on the unit circle of the complex plane, i.e., | A/| = 1 

and Xi 7 ^ X m for l ^ m. 

(2) : All eigenvalues are on the unit circle. There are repeated eigenvalues. But the geo¬ 

metric multiplicity for all eigenvalues are the same as the algebraic multiplicity, i.e., 
Mulg(Xi) = Mul a (Xi ) for all l. 

(3) : All eigenvalues are on the unit circle. There are repeated eigenvalues with algebraic 

multiplicity greater than geometric multiplicity, i.e., Mul g (Xi ) < Mul a (X{) for some l. 

(4) : There exits at least one eigenvalue not on the unit circle, i.e., |A;| ^ 1 for some l. 


Here, an eigenvalueA/ of M is a root of the characteristic polynomial Det(I — AM), and the 
algebraic multiplicity of an eigenvalue Mul a (Xi ) is the order of the root. The geometric mul¬ 
tiplicity of an eigenvalue Mul g (Xi ) is the number of independent eigenvectors corresponding 
to the eigenvalue. In general, Mul g (Xi ) < Mul a (Xi). According to the basic theory of linear 
algebra, Cases (1) and ( 2 ) are spectrally stable, and Cases (3) and (4) are spectrally unsta¬ 
ble. For Cases (1) and (2), we would like to know whether their spectral stability will be 
sustained under a small structural perturbation. Case (1) can also be shown to be struc¬ 
turally stable by considering the symplectic nature of M. As the structure of the system is 
perturbed, the eigenvalues will move accordingly. However, they cannot move off the unit 
circle due to a small structural perturbation on M for Case (1). This is because for every 
eigenvalue A off the circle, there exits another eigenvaule 1/A, which is on the opposite side 
of the unit circle. If one of the eigenvalues of Case (1) were allowed to move off the unit 
circle, then there would be more than 2 n eigenvalues. This forbidden situation is illustrated 
in Fig. [T] for n = 2. When there are repeated eigenvalues for Case (2), the constraints on the 
locations of the eigenvalues do not prohibit the eigenvalues moving off the unit circle due 


to structural perturbations, which is the so-called Krein collision 12 


m as illustrated 


m 


Fig. [2] for n — 2. Krein collisions preserve the symmetry of the eigenvalue distribution with 
respect to the real axis and the unit circle, and are the only possible pathways in parameter 
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FIG. 1. Eigenvalues are forbidden to move off the unit circle for Case 1. Illustrated here is the 
case of n = 2. 

space for a spectrally stable system to become spectrally unstable. When this happens, the 
system is structurally unstable. What is more interesting is that not all possibilities in Case 
( 2 ) are structurally unstable. Case ( 2 ) needs to be further divided into two sub-categories: 

(2.1) For all repeated eigenvalues, the corresponding eigenvectors have the same sig¬ 
natures. 

(2.2) There is at least one repeated eigenvalue whose eigenvectos have different signa¬ 
tures. 


Here, the signature of an eigenvector 0 of M is defined to be the sign of its self-product 
(il’, il’) = 0*fJ0. Here, J is the 2 n x 2 n unit symplectic matrix, i.e., 



The product between two vectors 0 and 0 in general is defined to be 


( 1 ) 


(0,0) = 0*fJ0, 


( 2 ) 


h Q, 


where 0* is the complex conjugate of 0 T . Krein 12], Gel’fand and Lidskii [131 ] . and Moser 
fl4 ] proved that Case ( 2 . 1 ) is structurally stable, and that Case ( 2 . 2 ) is structurally unstable. 
This is the celebrated Krein-Gel’fand-Lidskii-Moser theorem. 

Let’s use the example of a ID uncoupled lattice (n = 1) to demonstrate the process of a 
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FIG. 2. Eigenvalues are allowed to move off the unit circle for Case 2 when there are repeated 
eigenvalues. This is the so-called Krein collision. Illustrated in (a) and (b) are two possibilities for 
the case of n = 2. For each possibility, the eigenvalue distribution on the left is before the collision, 
and the one on the right is after the collision. 

Krein collision. For a ID uncoupled periodic lattice, the one-period transfer map M is 



M = Sq PS 0 , 


P = 


Sn = 


cos © sin < 


• Sill 0 COS ( 


( Wn 1 0 


y -w 0 w 0 
[ T dt 


Jo P (t) 

P (t) = w 2 (f) , 


(3) 

(4) 

(5) 

( 6 ) 
(7) 


where w(t) is a matched envelope function, w 0 = u>(0) and zb 0 = u>(0) are the initial 
conditions for w(t), and (f) is the one-period phase advance. The envelope function w(t) is 
determined by the envelope equation 

w + n q (t)w — w~ 3 . (8) 
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Since M is similar to P, and the eigenvalues and signatures are preserved by a similarity 
transformation [see Eq. (152]) ]. the spectral and structural stability properties of the lattice is 
completely determined by the phase advance matrix P. The eigenvalues, eigenvectors, and 
signatures are 


A + 

= cos (p + i sin (f ), 


= (1a) T > o-+ = -1, 

(9) 

A_ 

= cos (j) — i sin </>, 

7 \)_ 

= (1,-*) T , (T- = l. 

(10) 


Evidently, the system is spectrally and structurally stable when (j) ^ tot, which corresponds 
to Case (1). As the system parameters varies, the phase advance </> changes. A Krein 
collision occurs at </> = rm, where A+ = A_ and a + = — cm. This is Case (2.2), which is 
structurally unstable. Starting from a stable lattice with a small phase advance, we can 
increase the focusing strength and thus the phase advance gradually. The system is stable 
until the phase advance approaches tt. 


III. GENERALIZED COURANT-SNYDER THEORY 


To prepare for the investigation of the stability properties for general coupled lattice, we 
briefly summarize here the generalized Courant-Snyder theory, a thorough description of 


which can be found in Ref. 17]. 


The linear dynamics of a charged particle relative to the fiducial orbit are governed by a 
general time-dependent Hamiltonian 37] of the form 


H = ^z T Az , A = 


hi ( t ) R (t) 

R ( t) T m~ l ( t ) 


( 11 ) 


Here, z = (x,y,p x ,p y ) T are the transverse phase space coordinates, and n(t), R(t) and 
m~ 1 (t) are time-dependent 2x2 matrices. The matrices n(t), m~ l (t), and A are also 
symmetric. In this general Hamiltonian, the quadrupole component is in the diagonal terms 
of the hi (t) matrix. The off-diagonal terms of n (t) contain the skew-quadrupole and dipole 
components. The solcnoiclal component and the torsion of the fiducial orbit 38] are included 
in the R(t) matrix. The symplectic matrix specifying the map between Zq and 2 = M(t)z 0 
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IS 


fen — wRrr^j + ^mi? r + w (k — RmR T ^j — (w T wmw T ^j — 0 . (14) 


M{t) = S~ 1 P T S 0 , (12) 

( W ~ T 0 \ 

S=\ , (13) 

\ (wR — w)m w I 

where subscript “0” denotes initial conditions at t — 0, and w is a 2 x 2 envelope matrix 

function satisfying the matrix envelope equation 

d (dw 
dt 

In Eq. (jT2p . P G 5p(4) f|*50(4) = 1/(2) is a symplectic rotation, which is the generalized 
phase advance, determined by 

0 — n 

/i 0 

H = [wmw T ^j . (16) 

Alternatively and preferably, the transfer matrix M(t) can be expressed in terms of a 
symmetric envelope matrix u(t), which is defined to be the symmetric component of w(t) in 
its polar decomposition, 

= (17) 

j3(t) = w T (t)w(t). (18) 


P = -P 


(15) 


In terms of u, the transfer matrix is 

M{t) = S^P-'S; 


wO 

-1 


Su = 


u ‘ 0 

(uR — Du — u)m u 


— T~ 1 


D = L 


(■umu — umu ) + u(Rm — mR T )u 


(19) 

( 20 ) 
( 21 ) 


where P u G *S'p(4) p| *S'0(4) = U( 2) is a symplectic rotation determined by the differential 
equation 


P u 

ftU 


-Pu 


0 ft u 

ftu 0 


(umu T 'j 


-l 



( 22 ) 

(23) 







In Eq. (121 ft . L v l is the inverse of Lyapunov operator defined as 


L V {X) = LX + 1L, 


( 24 ) 


for a symmetric, positive-definite matrix V. A detailed discussion of the Lyapunov operator 
can be found in Ref. fl7]. The envelope matrix u(t) = \J (3{t) is determined from the 
differential equation for (3, 


P = 2e - fg - g T f - fill - h T (3 + 2m -1 /? -1 ™" 1 
e = uDu + u 2 — uD 2 u — uDu 
f = uDu + uu , 
g = (rh — Rm + , 

h = (k — RmR T — Rm — Rm)m ~ l , 


(25) 

(26) 

(27) 

(28) 

(29) 

(30) 


Equation (1251) is a second-order ordinary differential equation for /3, since every term on the 
right-hand-side is a function of ft and $. 

For every t, M(t ) is specified by two n x n symmetric matrices (3 and 0, and a U(n) 
matrix P u . The dimension of M(t) is thus (n 2 + n )/2 + (n 2 + n )/2 + n 2 = n(2n + 1), as 
expected for symplectic matrices. Another important advantage of using the symmetric 
envelope matrix u over the unsymmetric envelope matrix w is that Eq. (1191) enables the 
application of advanced techniques of pre-Iwasawa decomposition and normal form to find 
a matched solution for the j3 matrix without using the (inefficient) shooting method. 


IV. STABILITY ANALYSIS AND MATCHED LATTICE FUNCTIONS 

For practical applications of coupled lattices, it is desirable to design a coupled lattice be¬ 
longs to Cases (1) and (2.1), which are both spectrally and structurally stable. As mentioned 
previously, the parameter space satisfying this condition for a coupled lattice is larger than 
that for an uncoupled lattice. The generalized Courant-Snyder theory described in Sec. lIIIl 
provides an effective tool to study the stability properties of coupled lattices. One of the 
important result is that if a matched solution for (3 exists, then the stability property of a 
general coupled lattice is completely determined by the phase advance matrix P u . Using a 
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matched solution for (3, the one-period map is 

M(T) — S~q P u (T) t S u0 , (31) 

which indicates that M(T ) is similar to the inverse of the phase advance P U (T) T , and 
thus has the same eigenvalues and multiplicity as P u (T) t . Because P U (T) T is a rotation, 
its eigenvalues are on the unit circle. Now we show that the phase advance P U (T) also 
determines the structural stability of M(T ). For an eigenvector 0 of M(T), S0o0 is an 
eigenvector of P u (T) t , and the product between the two eignervectors defined in Eq. (j2J) is 
preserved by the similarity transformation, i.e., 


(Suo'ijj, S u0 (p) = ^*S^ 0 iJS U Q(j) = 0 *iJ(j) = (0,0) . 


(32) 


Therefore, P U (T) T and M(T ) have the same eigenvalues, signatures, and thus structural 
stability properties. 

As in the case of a ID uncoupled lattice, matched solutions for f3 are much more preferable 
than unmatched solutions, because the lattice functions can be completely determined by a 
matched (3 solution in one lattice period, i.e., the lattice functions are periodic in terms of 
the lattice period. On the other hand, if an unmatched (3 solution is used, the (3 function 
and other lattice functions have to be solved in the entire time domain of 0 < t < oo. Does 
a matched (3 solution always exist? The answer is negative. What are the conditions for the 
existence of a match f3 solution? We prove now that for Cases (1) and (2.1), Eq. (|251) admits 
matched solution for B. The proof utilizes the technique of normal form and pre-Iwasawa 


decomposition 


35 


36] for symplectic matrices. 


First, let’s invoke the established result that for Case (1), M can always be transformed 
into the following normal form with a symplectic matrix A, 


M = ANA ' 1 , 

( 


N = 


Ri 


Ri 


Rr 


(33) 


(34) 


/ 


Ri = 


COS 4>1 sin 0; 
^ — sin 0/ cos (pi 


(35) 
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Obviously, N G Sp(2n)f)SO(2n) = U{n ) is a symplectic rotation. This fact is proved 
from the existence of a complete set of 2 n orthonormal eigenvectors (0j, 0 _/), (l = 1, 2, n), 
satisfying 


(0£) 0m) ftImi 

( 0 —h 0 — m) &imi 

(0/,0_ m ) = 0. 


(36) 

(37) 

(38) 


Here, -0; and 0_; = 0; are a pair of eigenvectors corresponding to the eigenvalues A [ and A_/ = 
Xi, respectively. Equations and (1571) state that 0; and 0_; have different signatures. 
The normal form is actually explicitly constructed. The transfer matrix A is given as u 




(39) 


where and rft are real and imaginary components of the eigenvector 0;, i.e., 0; = + ipi- 

We now show that for Case (2.1), such a set of orthonormal bases exists as well. For 
a repeated eigenvector A with Mul g ( A) = Mul a ( A), the corresponding eigenvectors span 
a subspace M\ of R 2n . Because the signature never vanishes in M\, we can always select 
a set of orthonormal bases for M\ through a Gram-Schmidt process. The subspace M_a, 
a complex conjugate image of M\, has the same structure except that the signature has 
the opposite sign. Therefore, for both Case (1) and Case (2.1), the normal form given by 
Eq. ([55]) exists. 

Second, we apply the pre-Iwasawa decomposition to the symplectic matrix A. According 
to the theory of Iwasawa decomposition [351 . [36], a symplectic matrix G can always be 
uniquely factored as 

/ 


G = P 



(40) 


where P G Sp(2n) fj*S'0(2n) = U{n ) is a symplectic rotation, and Y and Q are symmetric. 
The statement is true as well if the decomposition is defined alternatively to be 


G = 



(41) 
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Let the unique pre-Iwasawa decomposition of A be 


Then, the transfer matrix is M 
and f3 such that 


A = P A S A , 


Y 0 


( 42 ) 


Sa = , • (43) 

\Q YY -'h 

= S^P^NPaSa- We choose the initial conditions for f3 


Suo = S A , 


(44) 


and the solution of [3 will give the same transfer map 


M = S- 1 P-\S A = Sa'P^NPaSa- (45) 

Thus 

S^P- 1 = S-^P^NPa. (46) 

The uniqueness of pre-Iwasawa decomposition requires that 


S-' = Sa 1 , (47) 

P-' = PANP a . (48) 

Equations (144)1 and (1471) prove that S~ l = S~q, i.e., the (3 solution is matched. Thus we have 
proven the proposition that for Cases (1) and (2.1), Eq. (144)) admits matched solutions. This 
proof is a constructive one, which can be used as an effective method to find a matched solu¬ 
tion for (3. The conventional method for Ending matched solutions is the shooting method, 
where one takes an initial condition for (3 and solves the ordinary differential equation once 
in one period. In general, the solution is not matched, i.e., (3(T) ^ /3(0) or $(T) ^ /3(0). 
The shooting method requires one to estimate a better initial condition based on the size of 
the mismatch, and solve the differential equation again. This iteration is carried out many 
times until a matched solution is found. The new method suggested by the above construc¬ 
tive proof of the existence of matched solution only requires solving Eq. (1251) once with an 
arbitrary initial condition. We can construct the one-period map M(T) using any matched 
or unmatched solution of Eq. (1441) . then the eigenvectors of M(T) can be calculated. When 
the set of bases satisfying Eqs. (144D - (155)) exists, the initial condition for a matched solution 
is uniquely given by Eq. (144)) . This new method applies to both ID uncoupled lattices and 
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coupled lattices in higher dimensions. Of course, this procedure fails when the set of bases 
satisfying Eqs. (1551) - (1551) does not exist. However, for the desirable Cases (1) and (2.1), such 
a set of bases exists. Another practical implication of the existence of matched f3 solution 
is that when an matched solution for (3 cannot be found, the lattice must be unstable. 


V. CONTINUOUSLY ROTATING QUADRUPOLE LATTICES 


As an illustrative application of the theoretical formalism developed in Secs. lIIIl and IIV1 
we investigate here the stability properties of a continuously rotating quadrupole lattice 
8 ]. The Hamiltonian of a charged particle in such a lattice is |3l|, [34[ 


H = -z t Az , A = 
2 



K(t) = KqO 


cos(27rt/T) sin(27rt/T) 
— sin(27rt/T) cos(27r t/T) 


(49) 

(50) 


We will smoothly vary K q oT 2 and observe the movement of the eigenvalues of the one-period 
map. For a given value of K q0 T 2 , we find a matched (3 solution using the procedure described 


above. The calculation is carried out using a code developed in Ref. [34]. The matched (3 
solutions for the cases of n q0 T 2 = 8 and n q oT 2 = 9 are plotted in Fig. [5] The eigenvalue 
distributions for different values of K q0 T 2 are plotted in Fig.U For n q oT 2 = 8 and K q0 T 2 = 9, 
it belongs to Case (1). The two eigenvalues on the left half circle have different signatures, 
and move towards one another as n q0 T 2 increases. At K q0 T 2 = 7r 2 , these two eigenvalues 
collide at A = —1. Since their signatures are different, this is an unstable Krein collision and 
thus structurally unstable. Right after the collision at K q oT 2 = 10, these two eigenvalues 
moves off the unit circle, and lead to an unstable lattice. This is the scenario depicted 
in Fig. [21(b). These analyses can be straightforwardly carried out for any coupled lattice, 


such as the N-rolling lattice 


31] and the Mobius accelerator [9]. We note that to calculate 


the eigenvalues displayed in Fig.Gfl it is not necessary to find the matched (3 solutions. Any 
solution of the [3 matrix over one lattice period can be used. Matched solutions are preferable 
when the (3 matrix is used as a lattice function for the machine. 
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FIG. 3. Matched solutions for the cases of (a) K q oT 2 = 8 and (b)«; g o T 2 = 9. 

VI. CONCLUSIONS AND FUTURE WORK 

We have studied in this paper the spectral and structural stability of charged particle 
dynamics in a coupled focusing lattice as a Hamiltonian system. The recently developed 
generalized Courant-Synder theory for coupled lattices has been applied. It is has been 
demonstrated that for coupled lattices that are spectrally and structurally stable, the matrix 
envelope equation must admit matched solutions. A new method is presented to determine 
a matched solution for the matrix envelope using the technique of normal form and pre- 
Iwasawa decomposition. If a matched solution exists, this method is able to determine the 
matched solution simply by solving the envelope equation once without using the (ineffi¬ 
cient) shooting procedure. As an example, stability properties of a continuously rotating 
quadrupole lattice are investigated. The Krein collision process for destabilization of the 
lattice has also been demonstrated. 


14 







1.5 r . . 1 -1.5 r . . 1 
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Re Re 


FIG. 4. Eigenvalues of continuously rotating quadrupole lattices. The two eigenvalues on the left 
half circle have different signatures in (a), and move towards one other as n q oT 2 increases in (b). 
The Krein collision occurs at n q oT 2 = n 2 in (c), after which these two eigenvalues on the left move 
off the unit circle and lead to an unstable lattice in (d). 


The application of coupled lattices to high-intensity charged particle beams 


39 


43] will 


be investigated in future studies. The theoretical framework and analytical tools devel¬ 
oped in the present study are also expected to be effective for the current theoretical and 
experimental investigations of emittance exchange technologies 


44 


50)]. 
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